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Abstract. Plasticity of two-dimensional discrete dislocation systems is studied. It is shown, that at some threshold stress level 
the response becomes stress-rate dependent. Below this stress level the stress-plastic strain relation exhibits power-law type 
behavior. In this regime the plastic strain rate induced by a constant external stress decays to zero as a power-law, which stems 
from the scaling of the dislocation velocity distribution. The scaling is cut-off at a time only dependent on the system size and 
the scaling exponent depends on the external stress and on the initial correlations present in the system. These results show, 
that the dislocation system is in a critical state everywhere we studied below the threshold stress. 
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INTRODUCTION 

Plastic deformation of crystalline materials is usually achieved by the motion of lattice dislocations. These linear 
crystal defects interact via long-range anisotropic stress-fields, and at low temperatures their motion is constrained to 
a glide plane. In a simplified model used in this paper only parallel edge dislocations are considered with a single glide 
plane. Although in crystals dislocations usually form complex three-dimensional networks, this model proved capable 
of reproducing many experimentally observed phenomena related to plasticity, like, e.g., strain avalanche statistics [1], 
Andrade-creep exponents [2, 3], and properties of X-ray profiles [4, 5]. 

This two-dimensional (2D) model was also found to exhibit a yielding transition, that is, below some yield stress Ty 
the plastic strain rate decays to zero, otherwise it tends to a constant value [2]. This fundamental observation raised the 
analogy with second-order phase transitions with Ty characterizing the critical point; for a review see [6]. This theory 
was elaborated by Laurson et al. by introducing a dynamical correlation length, that diverges as Ty is approached from 
above [7]. In addition, with a cellular automaton technique Zaiser et al. showed, that as the applied stress approaches 
Ty from below, both the total accumulated plastic strain and the cut-off of the strain avalanche distribution diverge 
[8, 6]. According to these results, the concept of a well-defined critical yield point seems to be established for this 2D 
model. 

In this paper it is shown, that a threshold stress Tth can be introduced with the plastic response showing power-law 
behavior below it, and strong dependence on the applied stress rate above it. It is argued, that below Tth the system 
behaves like one in a critical state, i.e., the characteristic cut-off times are diverging with increasing system size. This 
behavior was already confirmed in several different set-ups [9]. Here it is shown, that if the external stress is applied 
to the random/uncorrelated system, the exponents of the power-law relaxation change significantly, confirming their 
dependence on the initial conditions and external stress. According to these results, the concept of the existence of a 
non-equilibrium critical yield point for this 2D system is challenged, rather criticality for all relaxation states below 
Tth is suggested. 

THE DISLOCATION DYNAMICS MODEL 

A set of parallel straight edge dislocations with parallel slip planes is considered with a 2D representation on a plane 
perpendicular to the dislocation lines. For the motion of dislocations overdamped dynamics is assumed because of the 
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large acting friction forces. By introducing the notations r, — (x; ,y,) for the position of the rth dislocation, bi — Si{b, 0) 
for its Burgers vector (i, = ±1 is called its sign), the equation of motion of each dislocation takes the form [2] 



= 1; Jy^i 



(1) 



Here Tind('') = cos((j!))cos(2(j[>)r^ is the long-range shear stress field generated by an individual dislocation. Text is the 
external shear stress, and is the total number of the dislocations in the system. In the rest of this paper the different 
physical parameters are absorbed in the length-, time-, and stress-scales, as we measure them in the natural units of 



-0.5 



[pMGb ) , and Gbp , respectively, where p is the dislocation density, M is the dislocation mobility, and G 
is an elastic constant [10]. At the borders of the square-like simulation area periodic boundary conditions are applied. 



YIELD TESTS 

In order to study the plastic response of this system stress-controlled yield tests were performed [11]. The simulations 
were started from a random arrangement of an equal number of positive and negative sign dislocations (with a total 
number ofN= 128), which were then let to relax at zero external stress. After this the applied stress Text was gradually 
increased with constant rate. As seen in Fig. 1(a), individual simulations exhibit fluctuating stress-strain curves, with 
steps corresponding to sudden bursts of activity. One can, however, define an average stress-strain curve over an 
ensemble by assigning for every applied stress level Text the average of the plastic strain jpi values measured in the 
individual simulations [the result is the black thick line in Fig. 1(a)]. It was shown, that for small external stresses 
this curve is a power-law for at least two orders of magnitude [11]. At some threshold stress level Tih the power-law 
relation smoothly breaks down. Further analysis showed, that in the range of this Tth other characteristics, like the 
plastic strain rate, the fluctuation of the plastic strain, and the coefficient of the inverse cubic tail of the dislocation 
velocity distribution also behave similarly [11]. According to numerical fitting, an approximate value of Tth ~ 0.17 
was suggested [11]. It should be emphasized, however, that since the observed power-law breakdown is smooth Tth is 
not characterizing a single point, rather a transition regime. 

In Fig. 1 results about the stress rate dependence of the plastic response are reported. The averaging procedure for 
the stress-plastic strain curve shown in Fig. 1(a) and described above was repeated for simulations with a halved stress 
rate of text = 9 • 10^^ (in natural units). The result in Fig. 1(b) shows that the two average stress-strain curves overlap 
for small stresses, and split for larger stresses. As expected, at a given Text larger strains are observed for the smaller 
rate. The double logarithmic plot of Fig. 1(c) [with the same curves as in Fig. 1(b)] confirms the power-law type 
behavior below the threshold stress Tth. Figure 1(d) plots the difference of the two curves of Fig. 1(b) with switched 
axes, i.e., A7pi(Text) denotes the difference of the two average plastic strain values 7pi measured for the two strain rates 
at a given Text. (Note that in [11] the notation of Ay was used for a different quantity, the strain fluctuations.) It is 
seen, that Aypi also exhibits a transition at the threshold stress Tth, and the system is much more sensitive to the driving 
rate above Tth than below. In conclusion, the threshold stress introduced separates two distinctiy different regimes, 
presumably marking a yielding phenomenon. 



RELAXATION TESTS 

In this section the power-law regime below the threshold stress Tth is investigated. The following different simulation 
scenarios were suggested for this purpose: (i) relaxation of initially randomly distributed dislocations to an equilibrium 
state at zero external stress, (ii) response of the equilibrated system to an inserted fixed dislocation, and (iii) response 
of the equihbrated system to an external stress below Tth [9]. In all three cases relaxation to a steady state is observed 
with the following properties [9]; 

• The mean absolute velocity of the dislocations \v{t)\ in all scenarios and the plastic strain rate 7pi(f) in scenario 
(iii) decay to zero as power-law. This scaling regime is cut-off at a time ti only dependent on the system size as 
ti ~ v^, and not depending on the external stress. 

• In the scaling regime both the symmetric Pg and the antisymmetric part of the dislocation velocity distribution 
exhibit scaUng property as follows: 

/'s(v,0=?"/(?"v) and P,{v,t)=tT^g{tPv) (2) 
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FIGURE 1. Yield test with constant stress rate at a system size ofN = 128. (a) Stress-plastic strain curves for different realizations 
(thick color lines) and the average curve (thin black line), (b) The average stress-plastic strain curves for two different stress rates, 
(c) The same curves as in (b) now on a double logarithmic plot. The power-law regime for small stresses breaks down around the 
threshold stress Ttt,. (d) The difference of the average plastic strain values Aypi measured at a given stress for the two different stress 
rates (i.e., it is the difference of the curves in panel (b) after switching the axes). 



with appropriate / and g scaling functions. The power-law dependence of the mean absolute velocity and the 
plastic strain rate can be directly deduced from the above scaling of and P^, respectively: |v(f)| ^ f^" and 
Yj,i{t) t^-^P [9]. In the last expression 2j3 — 7 is called Andrade exponent. 
• The exponents are different in each case suggesting that they depend on the initial conditions and that, therefore, 
they are not determined only by the type of the interactions and dynamics of the system. 

Here the last point is elaborated, by repeating simulation scenario (iii) (response to a small external stress) with the 
stress applied not on an initially relaxed configuration, but on a completely random system. The applied external stress 
is Text = 0.17 ss Tth (at smaller stresses similar behavior is found). It is seen in Fig. 2(a) that the mean absolute velocity 
decays as |v(f)| ~ f^" with an exponent of a sa 0.72. The cut-off time behavior is identical to the one found before 
and described above. As in the other cases, the scaling of Ps is observed [Eq. (2)] in the scaling regime with the same 
a exponent [Fig. 2(b)]. Power-law decay characterizes 7pi(f ), too [Fig. 2(c)], also accompanied by the scaling of P^ as 
in Eq. (2) with j3 « 0.75 and 7« 0.5 [Fig. 2(d)]. The Andrade exponent of the strain rate is therefore 2j3 — 7« 1.0. 
Note, that in [10] with the same simulation set-up a stretched exponential form was obtained for 7pi(f), presumably 
because of too small system sizes. 

It is interesting to compare the observed exponents with the previous results. For a here 0.72 is found, while in 
scenario (i) (i.e., the same simulation with zero external stress) 0.85 was observed. So the scaling exponent is changed 
by only modifying the level of the external stress. On the other hand, when comparing the simulation of Fig. 2 with 
scenario (iii) (stress applied on a relaxed system) the only difference is in the initial configuration (relaxed with internal 
correlations or random). The exponents are again clearly different: j3 changes from 0.5 to 0.75 and 7 from to 0.5. 
The corresponding Andrade-exponent of 7pi(f ) is also shifted from around 0.6 to around 1 .0. The exponents, therefore, 
depend both on the level of the external stress and the statistical properties of the initial configuration. 



DISCUSSION 



It was shown, that below the threshold stress this 2D dislocation system exhibits slow relaxation with a cut-off time 
diverging with the system size. This behavior is typical for systems being in a critical point, it is, therefore, concluded 
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FIGURE 2. Relaxation from a random initial state with an applied stress of Text = 0.17. (a) Evolution of the mean absolute 
velocity for different system sizes A', (b) Scaling of the symmetric part of the velocity distribution Ps as Eq. (2) with a = 0.72 for 
the system size of N = 2048. (c) Evolution of the plastic strain rate fpi for different system sizes N. (d) Scaling of the antisymmetric 
part of the velocity distribution Pa as Eq. (2) with j5 = 0.75 and y = 0.5 for the system size of yV = 2048. 



that this system is always critical in this regime. Consequently, description of the yielding transition with a single 
critical point is challenged. In addition, the scaling exponents are sensitive to the initial conditions and to the external 
stress level. In particular, the Andrade-exponent 2j3 — 7 of the plastic strain-rate 7pi(f) is not universally 2/3 in this 
system, but can be as large as 1 depending on the properties of the starting configuration and the external stress. 

Further questions such as how this critical behavior is changed above the threshold stress and how Tth relates in fact 
to the yield stress Ty discussed in the Introduction still remain to be addressed. 
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